Influence of molecular rebinding on the reaction rate of complex formation.

We simulated Brownian diffusion and reaction-diffusion processes to study the influence of molecular rebinding on the reaction rates of bimolecular reactions. We found that the number of rebinding events, Nreb, is proportional to the target's size and inversely proportional to the diffusion coefficient D and simulation time-step Δt. We found the proportionality constant close to π-1/2. We confirmed that Nreb is defined as a ratio of the activation-limited rate constant ka to the diffusion-limited rate constant, kD. We provide the formula describing the reactivity coefficient κ, modelling the transient-native complex transition for the activation-controlled reaction rates. We show that κ is proportional to (D/Δt)1/2. Finally, we apply our rebinding-including reaction rate model to the real reactions of photoacid dissociation and protein association. Based on literature data for both types of reactions, we found the Δt time-scale. We show that for the photodissociation of a proton, the Δt is equal to 171 ± 18 fs and the average number of rebinding events is approximately equal to 40. For proteins, Δt is of the order of 100 ps with around 20 rebinding events. In both cases the timescale is similar to the timescale of fluctuation of the solvent molecules surrounding the reactants; vibrations and bending in the case of photoacid dissociation and diffusional motion for proteins.

The Hinge Region Strengthens the Nonspecific Interaction between Lac-Repressor and DNA: A Computer Simulation Study.

LacI is commonly used as a model to study the protein-DNA interaction and gene regulation. The headpiece of the lac-repressor (LacI) protein is an ideal system for investigation of nonspecific binding of the whole LacI protein to DNA. The hinge region of the headpiece has been known to play a key role in the specific binding of LacI to DNA, whereas its role in nonspecific binding process has not been elucidated. Here, we report the results of explicit solvent molecular dynamics simulation and continuum electrostatic calculations suggesting that the hinge region strengthens the nonspecific interaction, accounting for up to 50% of the micro-dissociation free energy of LacI from DNA. Consequently, the rate of microscopic dissociation of LacI from DNA is reduced by 2~3 orders of magnitude in the absence of the hinge region. We find the hinge region makes an important contribution to the electrostatic energy, the salt dependence of electrostatic energy, and the number of salt ions excluded from binding of the LacI-DNA complex.

Method for the analysis of contribution of sliding and hopping to a facilitated diffusion of DNA-binding protein: Application to in vivo data.

DNA-binding protein searches for its target, a specific site on DNA, by means of diffusion. The search process consists of many recurrent steps of one-dimensional diffusion (sliding) along the DNA chain and three-dimensional diffusion (hopping) after dissociation of a protein from the DNA chain. Here we propose a computational method that allows extracting the contribution of sliding and hopping to the search process in vivo from the measurements of the kinetics of the target search by the lac repressor in Escherichia coli [P. Hammar et al., Science 336, 1595 (2012)]. The method combines lattice Monte Carlo simulations with the Brownian excursion theory and includes explicitly steric constraints for hopping due to the helical structure of DNA. The simulation results including all experimental data reveal that the in vivo target search is dominated by sliding. The short-range hopping to the same base pair interrupts one-dimensional sliding while long-range hopping does not contribute significantly to the kinetics of the search of the target in vivo.

Fractal trace of earthworms.

We investigate a process of random walks of a point particle on a two-dimensional square lattice of size n×n with periodic boundary conditions. A fraction p≤20% of the lattice is occupied by holes (p represents macroporosity). A site not occupied by a hole is occupied by an obstacle. Upon a random step of the walker, a number of obstacles, M, can be pushed aside. The system approaches equilibrium in (nlnn)(2) steps. We determine the distribution of M pushed in a single move at equilibrium. The distribution F(M) is given by M(γ) where γ=-1.18 for p=0.1, decreasing to γ=-1.28 for p=0.01. Irrespective of the initial distribution of holes on the lattice, the final equilibrium distribution of holes forms a fractal with fractal dimension changing from a=1.56 for p=0.20 to a=1.42 for p=0.001 (for n=4,000). The trace of a random walker forms a distribution with expected fractal dimension 2.